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Abstract — This paper presents a scheme to design a tracking 
controller for a class of uncertain nonlinear systems using a 
robust feedback linearization approach. The scheme is com- 
posed of two steps. In the first step, a linearized uncertainty 
model for the corresponding uncertain nonlinear system is 
developed using a robust feedback linearization approach. In 
this step, the standard feedback linearization approach is used 
to linearize the nominal nonlinear dynamics of the uncer- 
tain nonlinear system. The remaining nonlinear uncertainties 
are then linearized at an arbitrary point using the mean 
value theorem. This approach gives a multi-input multi-output 
(MIMO) linear uncertain system model with a structured 
uncertainty representation. In the second step, a minimax linear 
quadratic regulation (LQR) controller is designed for MIMO 
linearized uncertain system model. In order to demonstrate the 
effectiveness of the proposed method, it is applied to a velocity 
and altitude tracking control problem for an air-breathing 
hypersonic flight vehicle. 

I. Introduction 

In this paper, a robust tracking control scheme is designed 
for a class of uncertain nonlinear systems. The design 
is composed of two steps. In the first step a linearized 
uncertainty model for the uncertain nonlinear system is 
developed using a robust feedback linearization approach. 
The feedback linearization approach has many applications 
in the process control and aerospace industries. Using this 
method, a large class of nonlinear systems can be made 
to exhibit linear input-output behavior using a nonlinear 
state feedback control law. Once the input-output map is 
linearized, any linear controller design method can be used 
to design a desired controller One of the limitations of the 
standard feedback linearization method is that the model 
of the system must be exactly known. In the presence of 
uncertainty in the system, exact feedback linearization is not 
possible and uncertain nonlinear terms remain in the system. 

In order to resolve the issue of uncertainty after canceling 
the nominal nonlinear terms using the feedback linearization 
method, several approaches have been considered in the 
literature [3], [3]-[5], [8], [9], [11], [19]. Most of these 
approaches use adaptive or related design procedures to 
estimate the uncertainty in the system. In these methods, 
mismatched uncertainties are decomposed into matched and 
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mismatched parts. These methods typically require the mis- 
matched parts not to exceed some maximum allowable bound 
[12]. The existing results that are based on mismatched 
uncertainties either do not guarantee stabihty or require some 
stringent conditions on the nonlinearities and uncertainties to 
guarantee the stability of the system [3], [5]. 

In this paper, we approach the uncertainty issue in a 
different way and represent the uncertain nonlinear system 
in an uncertain linearized form. We use a nominal feedback 
linearization method to cancel the nominal nonlinear terms, 
and use a generalized mean value theorem to linearize the 
nonlinear uncertain terms. In our previous work [15], [18], 
the uncertain nonlinear terms are linearized using a Taylor 
expansion at a steady state operating point by considering a 
structured representation of the uncertainties. This lineariza- 
tion approach approximates the actual nonlinear uncertainty 
by considering only the first order terms and neglecting all of 
the higher order terms. In [18], the uncertain nonlinear terms 
are linearized using a Taylor expansion but an unstructured 
representation of uncertainty is considered. In both of these 
methods, the linearized uncertainty model was obtained by 
ignoring higher order terms. In [17], we introduced the 
linearization of nonlinear terms using using a generalized 
mean value theorem [6], [13] approach. This method exactly 
linearizes the uncertain nonlinear terms at an arbitrary point 
and therefore, no higher order terms exist. In [17], the upper 
bound on the uncertainties is obtained by using unstructured 
uncertainty representations. The bound obtained using an 
unstructured uncertainty representation may be conservative 
which may degrade the performance of the closed loop 
system. In order to reduce conservatism and to obtain an 
uncertain linearized model with a structured uncertainty 
representation, a different approach for obtaining an upper 
bound is presented here. In contrast to [17], here we propose 
a minimax linear quadratic regulation (LQR) [14] controller 
which combines with a standard feedback linearization law 
and gives a stable closed loop system in the presence of 
uncertainty. Here, we assume that the uncertainty satisfies a 
certain integral quadratic constraint (IQC). 

The paper is organized as follows. Section presents a 
description of the considered class of uncertain nonlinear 
systems. Our approach to robust feedback linearization is 
given in Section HU] Derivation of the linearized uncertainty 
model and tracking controller for an air-breathing hypersonic 
flight vehicle (AHFV) along with simulation results are 
presented in Section |IV] The paper is concluded in Section 
|V]with some final remarks on the proposed scheme. 



II. System Definition 

Consider a multi-input multi-output (MIMO) uncertain 
nonlinear system with the same number of inputs and outputs 
as follows: 

m 

= f{x,p) + V gk{x,p)uk{t) + egip, x, u), 



k=l 



(1) 



\ v{t) 



yi{t) = Vi{x), i = 1,2,--- ,m 

where x{t) G R", u{t) = [ui uj^^ € 

[yi----ym\^ £ M™ and e 7^ 0. Furthermore, the system has 
norm bounded uncertain parameters lumped in the vector 
p G W^. Also, f{x,p), gi{x,p), and Ui{x,p) for i = 1, - ■ ■ ,m 
are assumed to be infinitely differentiable (or differentiable 
to a sufficiently large degree) functions of their arguments. 
The term g{p,x,u) in ([T]) is a nonlinear function which 
represents the couplings in the system. The full state vector 
X is assumed to be available for measurement. 

III. Feedback Linearization 

In this section, we first simplify the model ([T]i so that 
the term involving e vanishes. Here we assume that |e| is 
sufficiently small and hence indicates weak coupling. In gen- 
eral, e depends on the physical parameters of the system and 
may be available through measurement or known in advance. 
Instead of neglecting this coupling in the control design, we 
model eg{p,x,u) as an uncertainty function g{p,p,x) with 
certain bound, where, p denotes a new uncertainty parameter 
whose magnitude is bounded. The parameter p appears due 
to the removal of coupling terms which depend on the input 
u. Now we can write ([T]i as follows: 



x{t) = f{x,p) + ^gk{x,p)uk{t) 



fe=i 



(2) 



where f{x,p) — f{x,p) + g{p,p,x) is an infinitely dif- 
ferentiable function and p = [p p]'^ ■ Also, note that in 
equation which includes the new uncertain parameter 
p, we can write the system in terms of a new uncertainty 
vector p ^ pq + Ap G M^, where q — q + 1. Here, po is 
the vector of the nominal values of the parameter vector p 
and Ap is the vector of uncertainties in the corresponding 
parameters as follows: 



Po = [ Pio P20 ■ 
Ap = [ Api Ap2 



P{q-1)0 PqO \ , 
■ Ap(^_i) Apq 



We assume that a bound on \Aps\ is known for each s G 
{1, • • ■ ,q}. We also assume that the functions in the system 
(|2|l are differentiable. The standard feedback linearization 
method can be used on the nominal model (without un- 
certainties) by differentiating each individual element i/i of 
the output vector y a sufficient number of times until a 
term containing the control element u appears explicitly. The 
number of differentiations needed is equal to the relative 
degree of the system with respect to each output for 
i = 1,2, ■ ■ ■ TO. Note that a nonlinear system of the form 



^ with TO output channels has a vector relative degree 
r — [ri r2 ■ ■ ■ r^] [10]. We assume that the nonlinear 

m 

system (|2| has full relative degree; i.e. — where n 

is the order of the system. 

It is shown in [16] that in the presence of uncertainties ex- 
act cancellation of the nonlinearities is not possible because 
only an upper bound on the uncertainties is known: Indeed, 
we obtain 

m 

= fo{x,Po) +^gkQ(x,po)uk(t) 
fe=l 



Nominal part 



+ Af{x,p) 



k=l 



Agk{x,p)uk{t) 



(3) 



Unceitain part 

yi{t) = vi{x), « = 1,2,--- ,TO 

where. A/, and Ag^ are the uncertain parts of their respec- 
tive functions. After taking the Lie derivative of the regulated 
outputs a sufficient number of times, the system (|3]l can be 
written as follows: 



f*ix) +g^{x)u 



'rn 

L2fi'^i) + E L2-'[LAM)]uk 

^ k=l 



fe=l 



(4) 



where, 

f,{x) = [Lj^{u^)---LY^{um)f, 



fa 



and the Lie derivative of the functions with respect to 
the vector fields / and gk are given by 



_dui{x) 



f, L^jUi = Lf{L^j Lg^{ui) 



dui (x) 



-9k- 



dx ■" f ' f - " 

Note that in equation (|4|, we have deliberately lumped the 
uncertainties at the end of a chain of integrators. This is 
because the uncertainties in y}, - ■ ■ ,yi^^^ will be included 
in the diffeomorphism, which will be defined in the sequel. 
This definition of the diffeomorphism is in contrast to [16], 
where the uncertainties in yj, ■ ■ ■ ,yl'~^ are assumed to be 
zero; i.e. they satisfy a generalized matching condition [19]. 
The feedback control law 



u = -g*{x) ^f^{x)+g*{x) ^v. 



(5) 



partially linearizes the input-output map (|4|i in the presence 
of uncertainties as follows: 



y* = 



Vi 



(6) 



Nominal part 



Uncertainty part 



where AWp{x,u,p) = L2M) + E L2g,'[LfM)]uk, 

^ fe=i 

y* — [2/1^ •■••2/m"]"^' ^ — [vi----VmV is the new control 
input vector Furthermore, we define an uncertainty vector 
AW,; which represents the uncertainty in each derivative of 
the z"* regulated output as 



AW,{x,u,p) 



AWl{x,u,p) " 
AWi{x,u,p) 



AWp{x,u,p) 
and write for i — 1,2, - ■ ■ ,mas given below. 



(7) 



n 





^Wl{x,u,p) 
_ AWl^(x,u,p) _ 



(8) 



Let us define a nominal diffeomorphism similar to the one 
defined in [17] for each partially linearized system in ^ for 
i = 1, • • • , m as given below: 

Xi = Ti{x,pa) =[ fyi- yic y» - y»c yj ■■ y?"^ ]^ ■ 

(9) 

Using the diffeomorphism (|9|l and system ([8]l, we obtain the 
following: 



X = Ax + Bv + AW{x,v,p), 



(10) 



where x{t) = [xfW, • • • G M" (fi = n + m), 

v{t) = [vi V2 ■ ■ ■ v„i]'^ e M™ is the new control input vector, 
AW{x,v,p) = [AW[{.),AWn-), ••■ ,Al^,^(.)f is a 
transformed version of AW(x, u,p) using (|9]) and AM^i(.) = 

KW(.),u'r(.), ■•■ , ^«^^(-)F for * = 1,2, ... , m. 
Also, 
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In our previous work, these uncertainty terms are lin- 
earized at a steady state operating point and all the higher 
order terms in states, control and parameters are ignored 
in order to obtain a fully Hnearized form for ( fTOb . In this 
paper, we adopt a different approach to the linearization of 
the uncertain nonlinear terms in ( fTot . Here, we perform the 
Unearization of AW{x, v,p) using a generalized mean value 
theorem [6], [7] such that no higher order terms exist. 

Theorem 1: Let : M" ^ M be a differentiable 

mapping on M" with a Lipschitz continuous gradient Vwp'' . 
Then for given x x(0) in R", there exists a vector 

= X + i(x - x(0)) with t e [0, 1], such that 

4'\x) - 4'Hxm) = y4'Hc,).{x - xm- (n) 

Proof: For proof, see [7]. ■ 
We can apply Theorem [T] to the nonlinear uncertain part 
of ( [Tol l. Let us define a hyper rectangle 



» = { 



X 

V 



Xi<Xi< Xi 

Vi <Vt < Vi 



(12) 



where Xiy ™d denote the lower bounds and Xiy ™d 
Vi denote the upper bounds on the new states and inputs 
respectively. For this purpose, the gradient of (.) is found 
by differentiating it with respect to x v at an arbitrary 
operating point Cij — (x, ii, p) for i — 1,2,.-. ,to and 
i — 1,2,... ,ri where, w-^] G 58, and p e O. We 

assume w\^\xo,vo,po) — 0, x(0) = 0, and v{{)) = and 
write w'f\.) as follows: 



Jf^{x,v,p)=w'^\.i).[x' 
Then Aiy(.) can be written as 

AW{.) = $ 



X 
V 



(13) 



(14) 



where 



A. Linearized model with structured uncertainty representa- 
tion 

The equation ( fTOl ) can be written in a linearized form 
using (fT4l) . Note that the matrix $ is unknown. However, 
it is possible to write bounds on each term w'^p [cij) in $ 
individually and represent them in a structured form. For this 



purpose, we define each individual bound as follows: 

pi = ma.x{\\w'[^\cn)\\), 



(2) 

P2 = max(||w'i '{ci2)\\), 

Cl2 



Pri 

Pri+i 

Pri+r2-\ hi"„ 



= max(||«;'^'K,J||), 
= max(|hf)(c2i)||), 



(15) 



= max(||u;'|^'"''(c„ir„ 



Using the definition in ( fTSl l and (fl4] i. the model (fTOb can be 
rewritten as 

X(t) = (A + CiAiE-i + C-2.l^iKi + ■ ■ • + a-i A^i^ri + • ■ ■ 

+ Cf,l^f,k^)x(t) + (B + Ci AiGi + CaAaGa + ■ ■ ■ (16) 
+ C'riA^iG'ri H VCni^riGn)v{t), 

where Cfc for fc = 1, 2, • • • , n is a n x 1 vector whose fcth 
entry is one and the other entries are zeros, A'^ for fc = 
1,2,-- - ,nis a 1 xn vector whose fcth entry is pj. and the 

other entries are zero, Gi, • • • , Gri-i , Gri+i, • • • , Gn-i = 



for i = 1, 2, • 
below: 



• • ,m and G(.) is a 1 x m vector as defined 

Gr\ — [pri pri ' ' ' pri ] ? 
Gri+r2 — [pri+r2 Pri+r2 ' ' ' Pri+r2]i 



^ri+r2H \-rm ~ [pri+r2'\ h'^m Pri-\-r2-\ hr^, ' ' ' Pri'\-r2'\ \-r^n 

(17) 

and ||Afe|| < 1. Using the above definitions of variables, 
we will write the system in a general MIMO form as given 
below: 



■cit)^Ax{t) + Bv{t) + J2CkCk{ty, 



k=l 



(18) 



Zk{t) = kkx{t) + Gkv{t)- fc = l,2,-- - ,n, 



where x{t) e M" is the state; Cfc(t) = l^k[kkX + Giv] e M 
is the uncertainty input, Zk (t) G M is the uncertainty output, 
v{t) G M™ is the new control input vector and y{t) e R™ is 
the measured output vector 

Theorem 2: Consider the nonlinear uncertain system ([T]) 
with vector relative degree {ri,-- - ,rm} si x — x{Q). 
Suppose also that f{x{Q),p) — and v{x{Q)) = 0. There 
exist a feedback control law of the form Q and coordinate 
transformation (|9]l, defined locally around a;(0) transforming 
the nonlinear system into an equivalent linear controllable 
system ( fTSl ) with uncertainty norm bound p for ( fTsT i in a 
certain domain of attraction if 

1) L^fVjh ^ and L^gL^fh ^ for all < i < - 1, 



/ 



2) ri + r2 + • • • + r„i = n, 

3) the matrix j3 = [g*{x,po)\~^ exists, 

4) The uncertainty satisfies |A|<1 VxgS,wG 
*B, and p e 9. 



Proof: The proof directly follows from the form of the 
feedback control law (O which cancels all the nominal non- 
linearities and the linearize the remaining uncertain nonlinear 
terms using the generalized mean value theorem at [x w]^ S 
05, and p G 0. Since the generalized mean value theorem 
allows us to write any nonlinear function as an equivalent 
linear function, which will be tangent to the nonlinear 
function at some given points, we can linearize the remaining 
uncertain nonlinear terms using the generalized mean value 
theorem. Finally, it is straight forward to write the entire 
uncertain nonlinear system ([T]i in the linear controllable form 
(fTSI l by finding the maximum norm bound p in ( fTSI ) on the 
linearized uncertain terms in the region being considered, 
[x -5]^ € *8, and p e 6. ■ 

IV. AHFV Example 

A. Vehicle Model 

The nonlinear equations of motion of an AHFV used in 
this study are taken from the work of Bolender et al [2] and 
the description of the coefficients are taken from Sigthorsson 
and Serrani [20]. The equations of motion are given below: 



V = 



r cos a - D 



h = V sin 7, 



gsm7, 7 = 



L + T sin a g cos 7 



mV 



V 



0-7, Q = Myy/Iy 



(19) 



See [2], [20] for a full description of the variables in this 
model. The forces and moments in actual nonlinear model 
are intractable and do not give a closed representation of the 
relationship between control inputs and controlled outputs. 
In order to obtain tractable expressions for the purpose 
of control design, these forces and moments are replaced 
with curve-fitted approximations in [20] which leads to a 
curve-fitted model (CFM). The CFM has been derived by 
assuming a flat earth and unit vehicle depth and retains the 
relevant dynamical features of the actual model and also 
offers the advantage of being analytically tractable [20]. 
The approximations of the forces and moments are given 
as follows in [20]: 

L^g5'CL(Q,5e,5c,Ari,Ar2), (20) 

D^gS'CD(Q,5e,<5c,Ari,Ar2), (21) 

Myy ^ ztT + qScCM^a, Se, Sc, An, Ara), (22) 

T « q[(l>CTAa, An, Moo) + Cria, An, M^, Aa)], (23) 

~gCivJa,5e,5c,An,Ar2], i = 1,2,3. (24) 

The coefficients obtained from fitting the curves are given as 
follows. These coefficients are obtained by assuming states 
and inputs are bounded and only valid for the given range. 
Here, we remove the function arguments for the sake of 
brevity: 

CL = Cia + Cl^S, + Ci^S, + C^"i An + C^"^ Ara + Cl, 
Cm = C^ja + Ci}S, + d'rS, + Cf/^ An + C^"' Ara + CS„ 



Cn =Cf+^^'^ {a + An)' + C^o+^^'\a + An) 



+ CS''^Q<5e + C^"^5n 



Ct.4> — Ct,^oi + Cj,^j 

f-2 



Ct = C^-'Ad + C?a + + C^"i An + C?., 



Civ. = C;^, Q + C^=, 5e + C^=, + C^;i An + C^;^ Ara + C^r, . 

Here, Moo is the free-stream Mach number, and q is the 
dynamic pressure, which are defined as follows: 



9 = 



p(h)V^ 



V 

Mo' 



(25) 



Also, p{h) is the altitude dependent air-density and Mq is 
the speed of sound at a given altitude and temperature. The 
nonlinear equations of motion have five rigid body states; i.e., 
velocity V, altitude h, angle of attack a, flight path angle 7, 
and pitch rate Q. The CFM also has 6 vibrational modes and 
they are represented by generalized modal coordinates rji. 
There are four inputs and they are the diffuser-area-ratio Ad, 
throttle setting or fuel equivalence ratio (p, elevator deflection 
6e, and canard deflection 6c- In this example, tracking of 
velocity and altitude will be considered. 

B. Feedback linearization of the AHFV nonlinear model 

1 ) Simplification of the CFM: The CFM contains input 
coupling terms in the lift and drag coefficients. We simplify 
the CFM in a robust way as presented in Section |in] so 
that the simplified model approximates the real model and 
the input term vanishes in the low order derivatives during 
feedback linearization. In the simplification process, we will 
first remove the flexible states as they have stable dynamics. 
A canard is introduced in the AHFV model by Bolender 
and Doman [1] to cancel the elevator-lift coupling using an 

ideal interconnect gain fcec ~ 77 which relates the canard 

deflection to elevator deflection {5c = kcc5e)- In practice, 
an ideal interconnect gain is hard to achieve and thus exact 
cancellation of the lift-elevator coupling is not possible. In 
the simplified model, we assume that the interconnect gain 
is uncertain with a bound on its magnitude and it also 
satisfies an IQC [14]. The drag coefficient is also affected 
due to the presence of elevator and canard coupling terms in 
the corresponding expression. We also model this coupling 
as uncertainty. The simplified expressions for lift, moment, 
drag, and thrust coefficients now can be written as follows: 



Cr 



Cm = Cfja + [C^J - C'mC— §^)]5e + C%j, 



(26) 



f-2 



■Cl 



Ct = Cf^'^Ai + C^a + C' 



where AC; is the uncertainty in the lift coefficient (7l(q!) 
due to the uncertain interconnect gain and AC^ is the uncer- 
tainty in the drag coefficient C^i (a) due to the input coupling 
terms. Furthermore, in order to obtain full relative degree for 
the purpose of feedback linearization, we dynamically extend 
the system by introducing second order actuator dynamics 
into the fuel equivalence ratio input as follows: 



(27) 



After this extension we have two more states 4> and 0, and 
thus the sum of the vector relative degree is equal to the 
order of the system n; i.e. n — 7 and thus satisfying one 



TABLE I: States, Inputs and Selected Physical Parameters at 
the Trim Condition 



Vehicle Velocity V 


7702.08 ft/sec 


Altitude h 


85000 ft 


Fuel-to-air ratio 


0.4388 


Pitch Rate Q 





Angle of Attack a 


—0.0134 rad 


Flight Path Angle 7 





Elevator Deflection <5e 


2.005 deg 


Canard Deflection 5c 


2.79 deg 


Diffuser Area ratio A^j 


1 


Reference Area S 


17 sq-ft.ft-^ 


Mean Aerodynamic Chord c 


17 ft 


Air Density p 


6.6 X 10--^ slugs/ft^ 


Mass with 50% fuel level m 


147.9 slug, ft-l 


Moment of Inertia lyy 


5.0 X 10^ slugs/ft2/(rad . ft) 



of the conditions for feedback linearization [10]. We use 
Theorem |2] to linearize the AHFV dynamics. The outputs 
to be regulated are selected as velocity V and altitude h 
using two inputs, elevator deflection Sc, and fuel equivalence 
ratio (pc- Note that the canard deflection Sc is a function 
of the elevator deflection Sc and they are related via an 
interconnect gain. Also, we fix the diffuser area ratio Ad to 
unity. The new simplified model consists of seven rigid states 



x=[Vhja4'4>Q 
by a general form as follows: 



and can be represented 



(28) 

where the control vector u and output vector y are defined 

as 

u = [ui,u2f = [SeAcf ,y = [yi,y2f = [v,hY' . 

The following set of uncertain parameters are considered 
for the development of a linearized uncertainty model: 



[CI 
ACm 



c" 



ACi ACd ACt 



AG 



T,4>\ 



(29) 



We assume that p G O, where 8 = {p G | 0.9po ^ Pi ^ 
l.lpo for i ~ I, ' - :9}. In order to get the linearized 
uncertain model for the uncertain nonlinear AHFV model, 
the output V and the output h are differentiated three times 
and four times respectively using the Lie derivative: 



V 



where 



= ,f.(a;,Po) + sr.(K,Po) 



Se 

4>c 



+ 



A V {x, u,p) 
Ah'^{x, u,p) 



Nominal nonlinear part 



Uncertain nonlinear part 



(30) 



L,,{L).h) 



and A V {x,u,po,Ap) and Ah'^{x,u,po, Ap) are the 
uncertainties in V and h'^ respectively. The application of 
the control law (|5]l yields the following: 



V 



VI 
V2 

Nominal linear part 



+ 



A V {x, u,p) 
Ah'^{x, u,p) 

Uncertain nonlinear part 



(31) 



- V ' 









V 









V 




VI 




h 







+ 


h 









h 









. . 




V2 





Also, by using the fact that there are no uncertainty terms 
in V, and h, we can write linearized input-output map for 
the AHFV model using (HJ as follows: 


AV 

A V {x,u,p) 

. (32) 

Ah 

A h {x,u,p) 
A/i^(a;, u, p) 

Let us define a diffeomorphism for each system as in (|9|l 
which maps the new vectors ^ and -q respectively to the 
original vector x as follows: 

C = ri(a;,po,V;), rj = T2{x,po,h,), (33) 

where 

Ti {x, p, Vc) = [ /J {V{t) - V,)dT V-Vc V V ]^ , 

T2{x,p,h,)=\^ p^{h{T)~K)dT h~K h h h^, 

and Vc and he are the desired command values for the veloc- 
ity and altitude respectively. We write each diffeomorphism 
as follows: 

X = T{x,po,Vc,hc), (34) 
where X = [ §2 ^4 Vi V3 Vt Vs and 

' n{x,p,Vc: ' 

T2{x,p,hc] ^ 

the nominal part of ( |32] | into the new states using the 
transformation ( l34b and linearize the uncertainty parts of OTI ) 
using the generalized mean value theorem as follows: 



T{x,p,Vc,hc 



Now, we can transform 



Kr = [ pr p7 P7 P7 PT P7 P7 P7 P7 \ , 

ks = [ PS pa PS p8 PS PS PS PS PS ] , 

^9 = [ P9 P9 P9 P9 P9 P9 P9 P9 P9 ] , 

Ki = K2 = ks = Ks = [ 0]; 

G4= [ P4, P4 ] , 
Gg = [ p9 P9 ] , 

Gi=G2=Gs = G5 = G(i = G7 = Gs=[ ] . 

Using the above definition of the variables, we can write 
the system in the general MIMO form given in ( fTSl l, where 
x{t) e MP is the state, (k{t) = Ak[KkX + Giv] e R is the 
uncertainty input, Zk{t) G M is the uncertainty output, and 
v{t) G is the new control input vector. 

C. Minimax LQR Control Design 

The linearized model ([18 1 corresponding to the AHFV 
uncertain nonlinear model (|l"9]l allows for the design of 
a minimax LQR controller for the velocity and altitude 
reference tracking problem. The method of designing a 
minimax LQR controller is given in [14]. Here, we follow 
the same method and proposed a minimax LQR controller 
for the linearized system ( flSl ). We assume the uncertainty in 
the system ( fTSl ) satisfies following IQC and the original state 
vector X is available for measurement. 



(II ^At) 11-11 iAt) \\ )dt > -X^(0)D.X(0), (37) 



where Dj > for each j = 1 , ■ • • , to is a given positive 
definite matrix. The cost function selected is as given below: 
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[xitfQxit)+vitfRvit)]dt, 



(38) 



where Q = > and R = > are the state and 
control weighting matrices respectively. A minimax optimal 
controller can be designed by solving a game type Riccati 
equation 



(39) 



In this section, we write the e quatio n ( 1351 ) in a structured 
form as presented in Subsection IIII-AI Using ( fTSl ), ( fT6l l, and 
( [TtI i we can write ( 1351 ) as given below: 



where 



Xit) = {A + CiAi A'l + C2A2A'2 + • ■ • + CgAgKg)x(t) 
+ {B + Ci AiGi + C'2A2G2 + ■ ■ ■ + CgAgGsMt) 
where 
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The weighting matrices Q and R, and parameters t^, for 
A; = 1, 2, • • • , 9 are selected such that they give the minimum 
bound 



min[x^(0)X,x(0) + ;^r,x^(0)73,x(0)], 



(40) 



on the cost function ( 1381 ). The minimax LQR control law 
can be obtained by solving the ARE ( |39] | for given values 
of the parameters as given below: 



V{t) = -GrX{t), 



(41) 



where 



Gr = E-'^[B^Xr + G'^K] 



is the controller gain matrix. The parameters Q and R 
are selected intuitively so that required performance can 
be obtained and for k = 1, 2, • ■ • ,9 correspond to the 
minimum bound on (l38T l. These parameters are given as 
follows: 



Q = diag 
R 



3.0 





1000, 500, 500, 100, 0.001, 100, 100, 500, 500 ] , 
T3 = 547.9, T4 = 8.0, TT = 4935.7, 
rg = 4935.3, rg = 3768.0. 





3.0 



D. Simulation Results 

The closed loop nonlinear AHFV system with the minimax 
LQR controller i4T[ is simulated using different sizes and 
combinations of uncertainty. For the sake of brevity, here 
we evaluate the performance of the proposed controller by 
using step input commands for the following three cases: 

1 . Uncertain parameters equal to their nominal values, with 
no uncertainty. 

2. Uncertain parameters 20% lower than their nominal 
values. 

3. Uncertain parameters 20% larger than their nominal 
values. 

The responses of the closed loop system given in Fig. [T] - 
Fig. [8] show that the minimax LQR controller along with the 
feedback Unearization law gives satisfactory performance. 



Fig. 2: Control input responses corresponding to a step 
change in velocity reference. 
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Fig. 3: Responses of a and 7 for a step change in velocity 
reference. 
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Fig. 1: System response to a step change in velocity refer- 
ence. 
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Fig. 4: Flexible states during velocity reference tracking. 



V. Conclusion 

In this paper, a robust nonlinear tracking control scheme 
for a class of uncertain nonlinear systems has been proposed. 
The proposed method uses a robust feedback linearization 
approach and the generalized mean value theorem to obtain 
an uncertain linear model for the corresponding uncertain 
nonlinear system. The scheme allows for a structured un- 
certainty representation. In order to demonstrate the appli- 
cability of the proposed method to a real world problem, 
the method is applied to a tracking control problem for an 
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Fig. 5: System response to a step change in altitude reference. 



20 30 



50 60 



50 60 



Fig. 6: Control input responses corresponding to a step 
change in altitude reference. 
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Fig. 7: Responses of a and 7 for a step change in altitude 
reference. 



air-breathing hypersonic flight vehicle. Simulation results for 
step changes in the velocity and altitude reference commands 
show that the proposed scheme works very well in this 
example and the tracking of velocity and altitude is achieved 
effectively even in the presence of uncertainties. 
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Fig. 8: Flexible states during altitude reference tracking. 
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